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Abstract 

We prove that the first passage time density p{t) for an Ornstein-Uhlenbeck process X{t) obeying 
dX = —PX dt + a dW to reach a fixed threshold d from a suprathreshold initial condition xo > > has 

a a lower bound of the form p(t) > k exp [— pe®''*] for positive constants k and p for times t exceeding some 
positive value u. We obtain explicit expressions for k,p and u in terms of /3, a, xo and 6, and discuss 
application of the results to the synchronization of periodically forced stochastic leaky integrate-and-fire 
I model neurons. 

> 

<^ 1 Introduction 
cn 

The distribution of the first passage time t of an Ornstein-Uhlenbeck process to a fixed boundary is of 
I interest in a variety of fields including mathematical neuroscience, where it describes the inter-spike interval 

• • distribution of a leaky integrate and fire (LIF) model neuron driven by a combination of steady and fluc- 

. tuating current sources [I] . Unlike the density of the first passage time to a fixed boundary for Brownian 

motion, the FPT density for the OU process is not known in closed form. Explicit expressions are known 
for the moments of the first passage time; in a neuroscience context the mean and variance are related 
respectively to the rate and variability of neuronal discharge [HI [28] . Conditions guaranteeing the existence 



and smoothness of the FPT density to constant and moving boundaries have been studied [HI [H] . 



Intuitively it is perhaps "obvious" that the FPT density for an OU process is never identically zero, because no 
matter how removed a time t is from the typical time of first passage, there must exist a set of trajectories 
of positive measure that make a sufficiently large excursion from the mean behavior that arrival at the 
boundary is delayed until t or later. It is the purpose of this article to augment this intuitive argument 
by providing an explicit positive lower bound valid for the tail of the distribution. We restrict attention 
to the lower bound problem for the "suprathreshold" case, i.e. the case in which the initial value xq and 
the asymptotic mean of the process are on opposite sides of the threshold value. In the neural context 
this case corresponds heuristically to that of an LIF model neuron driven by a sufficiently strong injected 
current that it will reach its firing threshold in finite time, i.e. firing is not dependent on the presence of 
fluctuations. This case is therefore fundamentally distinct from that of "stochastic resonance" phenomena, 
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which occur in a subthreshold setting. For the subthreshold case the initial and asymptotic value of the OU 
process lie on the same side of the capture boundary. In this case asymptotic results have been proven that 
characterize the tail of the distribution at large times [HITS]. In contrast, the author is aware of no results 
providing strictly positive lower bounds for the FPT density in the suprathreshold case. Such bounds would 
be of interest in the study of synchronization of neural activity via periodic stimulation, for reasons to be 
discussed below. Our main results are summarized as Theorem [T] and the following Remark. 

Theorem 1. Positive Lower Bound Theorem. Let {X[t), t>Q\ he an Ornstein-Uhlenbeck process satisfying 
the Ito stochastic differential equation 

dX ^ -(iXdt + adW (1) 

with initial condition 

X(0) = xo, (2) 

where P, a and xq are positive constants and W{t) is a standard Wiener process. Let 6 be a fixed threshold 
satisfying xq > 9 > 0. Then there are positive constants k,p and u such that the first passage time density 
of the process to the threshold, px(t), satisfies 

PxW >fcexp[-pe6^*] (3) 

provided t > u. 

Remark 2. Given particular values of the constants (3,a,xo and 9, the inequality Q) is satisfied for the 
following values of k,p and u: 



(f-l) (4) 
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where V denotes the maximum operator. 

The proof is obtained from a geometric construction, elaborated in Sections [2][3j that exploits the well known 
change of variables by which an OU process may be written in terms of a standard Brownian Motion |25) . 
Sato [H] used the same transformation to obtain analytic results on the FPT distribution of the OU process 
to a constant boundary, in the subthreshold case. While our results use the same change of variables, one 
cannot obtain the explicit lower bound sought directly from Sato's asymptotic results. 

The plan of the paper is as follows: in Section[2]we introduce the needed change of coordinates and establish 
a lower bound for the first passage time of a standard Brownian motion to a piecewise linear approximation 
of a square root boundary. In Section [3] we use these results to establish Theorem [T] In Section |4] we discuss 
potential applications of our result to synchronization of periodically forced leaky integrate and fire model 
neurons in the presence of additive noise, as well as related asymptotic results in the existing literature. 



2 Bounding the FPT distribution to a square root boundary for 
Brownian motion 



2.1 Coordinate Transformation 



Let X{t) be an OU process described by Equations 
i.e. B{t) satisfies 



Tp), and let B{t) denote a standard Brownian motion. 
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dB = dW (7) 
5(0) = 0, (8) 

from which E[B] — and Y[B] — t. We wiU identify the sample space ft for Brownian trajectories with 
C[0,cx)). 

For t > we have the representation of X in terms of B: 

X{t) = xoe-^* + '^B (e^'^* - l) . (9) 



The problem of finding the first passage time distribution for an OU process with a fixed boundary is 
equivalent to the problem of finding the FPT distribution for a standard Brownian motion to a moving 
boundary with square root time dependence. Given the parameter (3 for the OUP, Equation ([l]), define an 
exponentially rescaled time 

s = e^^* - l,t > 0. (10) 



Equivalently, for s > 0, we have t = (^jp^ log(s + 1). 



Lemma 3. Let X and B be as defined in Equations ([TH^, with xq > 9 > 0, and let s be as defined in 
Equation (lO^ . Fix a particular Wiener process trajectory w G O. Then X{t,uj) > 9 for < t < t if and 
only if 

j2B 

B{s, u)> {9^/7+l - .To) 

for allO<s< e^^^ - 1. 

Proof of Lemma From Equation ^ , 

X{t) > 9 for < i < T (11) 

^ xoe-'^* + ^^B {e^P* -l)>9 for < t < r (12) 

^ B {e'P' - 1) > {9 ~ xoe-^*) e'^* for < i < r (13) 

<^ B{s) > {e^'9 - xo) for < s < e^'^^ - 1 (14) 

B{s) > ^ {9y/7TT ~ xo) for < s < e^^^ - 1. (15) 

This concludes the proof of Lemma [S] □ 

Note that when s = t — 0, we have B{0) = > ■\J~^{9 — xq), consistent with the assumption that xq > 9. 
Let t' be the time of first passage of X(t) to 9. Let s' be the time of first passage of B{s) to the boundary 



bis) = ^ {9V^^ - xo) . (16) 
The following corollary to Lemma [3] is immediate. 
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Corollary 4. Let the processes X, B be defined as in Lemma\3\ Let b be the boundary function given by 



Equation (16) and let s be the rescaled time given by Equation llw. Then as dt ~> the first passage time 
densities of X and B satisfy 

Pr{t' e[t,t + dt)) = Pr(s' e [s, s + ds)) 

with 

ds = dt = 2/?e^''* dt = 2f3{s + 1) dt. 

Consequently px{t), the first passage time density for X, can be obtained from the first passage time density 
for B, pb{s), as 

pxit) - 2/3(s + l)pBis) - 2(3e^^'pB (e'^* - l) . (17) 
An explicit positive lower bound for pB will then provide a positive lower bound for px- 



2.2 Piecewise Linear Approximation for b{s) 

In order to obtain a positive lower bound for pB '^a will need the slope 6'(s) of the boundary function b{s): 



^ ' \ <J J 2V^TT V 2 ^/^TT 

In particular, we have b'{0) = 9 /cr > and for all s > 0, 6'(s) > 0. 

Given s' > and a small time ds > let fip C be the set of all trajectories in the sample space fl — C[0, oo) 
with first passage time s' to the boundary b{s) in the small interval s' € [s, s + ds). We will decompose the 
set Oq into two subsets and obtain a positive lower bound for the measure of one of them. First we construct 
a piecewise linear approximation of the boundary function b (see Figure [T]) . 

Let Li be the half line {(s, ai + sbi)\s > 0}, where 

ai = b{0) = {0 -XQ)y/2/3/a <0 (19) 
6i = &'(0) = (6'/2)v/2/3/ct > 0. (20) 



The line Li is tangent to the boundary function b(s) at s = 0. The boundary function has negative second 
derivative for all s > 0, so we have ai + sbi > b{s) for all s > 0. Consequently every continuous trajectory 
B{s) that begins at the initial condition -B(O) = and meets the boundary b{s) at some time s' > must 
first meet the line Li at some time s" < s' . 

Let s' > be given. We define a second boundary line L2 as the horizontal line with height 



/2P 

= "[^ ) = - 

Lines Li and L2 intersect at the point 



a2 ^ b{s') = ^ i^eV^^ - xoj (21) 



,^^^l^=2(VVTl-l). (22) 



bi 

Note s* < s' provided s' > 0. The duration A of the interval I — [s*, s'] is 



A = s' - s* = s' + 2(1 - Vs' + 1) (23) 



4 




1 2 3 4 5 

Time s 



Figure 1: Geometric construction of the piecewise linear boundary approximation, and sample trajectories. 
Here the boundary b{s) is the arc starting from 6(0) « —2.8, and is given by Equation (16) with a = .5, (5 = 
1, xo = 2 and 9 = 1. The line Li begins at the same location and rises above the arc with slope bi = Line 
L2 is horizontal, intersecting line L-y at = 2\/h — 1 2.47 and the boundary 6(s) at s' = 4. Approximate 
Brownian trajectories were generated using Matlab beginning at -B(O) = and evolving with noise factor 
o = .5. Of the twenty-five trajectories shown, three (thicker lines, shown in red and cyan) cross L\ for the 
first time at s" > s^,. One of these trajectories (thickest line, shown in cyan) then crosses L2 between in the 
interval [s' , s' + ds); here ds = 0.2 for illustration. Because L2 lies below b{s) for s > s' , the cyan trajectory 
must also have its first passage to b{s) within the interval [s',s' + ds). The cyan trajectory is a typical 
element of fl2. 
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and satisfies A + 2\/A = s' . Figurcjljillustrates the geometry of 6(s), Li and L2- We note for future reference 
that when s' — 8, s^, — A = 4, and that s* is a strictly increasing function of s' for s' > 0. 

Given ds > 0, we have defined Hq to be the collection of all Wiener process trajectories oj E fl = C[0, oo) 
such that the process X{s, uj) makes its first passage from X{0) = to b{s) in the interval [s', s' + ds). We 
further define fli to be the subset of flo consisting of those trajectories in Qq that make their first passage 
from X{0) = to the line Li at any time s" > s*. Figure [I] shows an example of such a trajectory highlighted 
in cyan (light blue). Starting from X(0) = the cyan trajectory remains above the line Li until passing the 
intersection of Li and L2 at s = s*, after which it crosses b{s) for the first time in the interval [s', s' + ds). 
Of the trajectories in fii, some - like the trajectory shown in cyan - also cross the line L2 somewhere in the 
interval [s', s' + ds). Denote the set of trajectories in fli with this property by ^2- 

Clearly, C ili C ^Iq. Because 6(s) is continuous and strictly increasing, every trajectory that remains 
above line Li until at least s = s* and then intersects line L2 for the first time in the interval [s',s' + ds) 
must also make its first passage to b{s) in the interval [s',s' + ds). Hence if the density of trajectories in 
^2 crossing L2 in [s', s' + ds) is strictly positive, so is the first passage time density of all trajectories in flo, 
which is the PPT density of interest. We next derive an explicit positive lower bound on the PPT density 
of trajectories in il2, hence also a lower bound for the PPT density of trajectories in ^Iq. 



2.3 Strictly positive lower bound for pB 



The first passage time density for a standard Brownian motion from B{0) = to a boundary moving with 
constant velocity is well known {e.g. [3], p. 221; [H]). The density of first passage from B{0) = to Li at 
time s is 

/ N f iai+bis)^\ 
9o^(^)= (^^^ -P [ ^ ) (24) 

where ai and bi are given by Equations ( T9]|20 ). Similarly, the density of first passage from a starting 
position fli + bis at time s to the line L2 at time s' > s is given by 

, n ^ I(qi +bis)-a2\ ^ ( ((ai - 02)^ ^ , , 

Because of the almost sure continuity of Brownian trajectories, we may write the first passage time density 
172(5') for trajectories in rt2 from initial position B{0) = to line L2 at time s' as a convolution of goi with 
gi2. Using li{s) to denote (oi + bis), we write 

92{s') = / goiis)gi2{s'\s)ds (26) 



Js, ((s'-s)s)3/2 2s 2(s'-s) 

where the upper end point (s') and lower end point (s*) for the integration both depend on s', as does the 
height of L2, namely 02 = b{s'). The integrand in Equation ( |27[ ) approaches zero at both endpoints, but is 
stricly positive in the interior of the interval / = {s|s* < s < s'}. To obtain a positive lower bound for the 
integral we rewrite the integrand as a product of three factors: 

92{s') = ^ / F,{s,s')F2{s,s')Fs{s,s')ds (28) 
271" Jsei 

Fi{s,s') = ai - a2{s') + bis (29) 
F2is,s') - (s{s' - s))-'/' (30) 
f is' ~ s)ll{s) + sjhjs) - a2)' \ f-/3 Qi{s,s') \ 
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where we introduce the notation 



Qi{s,s') = {s'll{s)^2sk{s)a2 + sal) {2cr^/^) (32) 

= s^9{Axo + e{s' - Ay/l + s')) + 4s'(a;o - Of (33) 
-4s(a;2 + (s' - 2)xoe + (-1 - 2s' + 2^/l^)e'^) 

Q2{s,s') = s(s'-s). (34) 

By construction, Fi(s*,s') = for ah s' > (recaUing that depends on s') and Fi(s,s') > for s > s*. 
On the other hand, i^3(s, s') — as s — )■ s'. Thus the integrand approaches zero at both ends of the interval 
of integration. In order to obtain a strictly positive lower bound on the integral we restrict integration to a 
subinterval of [s^, s']. Let rj, v be chosen in the interval ^ <t] <v < \ and set 

S\ = s^, + 77A (35) 
S2 = s* + (36) 

where as above A = (s' — s^) is the duration of the integration interval. Then < si < S2 < s'. Let /' 
denote the interval /' = [si,S2] C /. Each factor Fi,_F2,^3 is nonnegative on the interval /. Consequently 
we have the inequality 

ff2(s')>^/ F,{s,s')F2{s,s')F,{s,s')ds. (37) 
We consider each factor in turn. 



Factor Fi{s, s') is strictly increasing in s for all s, s' so for all s G /' we have 



Fi(s,s')>i^i(si,s') = r,^A. (38) 



The quadratic Q2(s, s'), which also appears in factor F2{s, s') = Q2is, s')^^^^ , has a global maximum (in s) 
of (s')^/4 at s = s'/2. In the range s € [si, S2], Q2(s, s') is strictly positive. Consequently for all s e /' we 
have 

F2is,s')>F2{s'/2,.s')^^. (39) 



From Equation (27) it is clear that both Qi{s, s') and Q2{s, s') are strictly positive for s G /'. The factors 
/3 and a are also positive, therefore for s e /', 

^3(s,sO>exp(-A "---^'gf (40) 
\ 4a-"= mm// (^2(5, s ) / 

As established in the next Lemma, Qi is just a quadratic function of s with second derivative that is 
guaranteed to be positive, provided s' > 8. 



Lemma 5. // s' > 8 and xq > 6 > then the quadratic function Qi(s,s') defined by Equation (32) has 
positive second derivative with respect to s. 



Proof. The second derivative of Qi with respect to s is constant in s, d^Qi/ds^ — 26{Axq + 6{s' — A^/l + s')). 
Suppose s' > 8. Then (s')^ > 8s', so (4 + s')^ > 16 + 16s'. Taking square roots, we obtain 4 + s' > 4\/l + s', 
which means that 1 > \/l + s' — s'/4, and since > by assumption, 9 > 6{\/l + s' — s'/4). By the 
assumption on xq, Xq > > 9{^/l + s' — s'/4). Consequently xq + 0{s' /A — ^/l + s') > 0, from which the 
conclusion follows. □ 
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Because Qi{s,s') has constant positive second derivative with respect to s, it achieves its maximum on /' 
either at si or at S2- The quadratic Q2{s,s') has constant negative second derivative with respect to s and 
achieves its minimum at one of the end points as well. Therefore we may conclude that[^ 



Fsis) > exp - 



4a2 



Qi(si) VQi(s2) 

Q2(si) AQ2(S2) 



(41) 



where V denotes the greater and A the lesser of the two arguments. 



Therefore in terms of the quadratics Qi, Q2, and the given constants, we have the following lower bound for 
the first passage time density at time s' 



(WW 

2a 



A 



4cr2 



Q2{si) A Q2{S2) 



(42) 



provided s' > 8. 

Note that for s' > 8, A(s') > |(s' - 2) > 4. Since |ai| = {xa - 9)y/2(3/a > we may write 
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g2{s') > {'y-vr—^_ (n^'?) 1 (^jexp^ 



27rcr 

64p9 f xo-d \ ( fi 



9 a 



Qi(si) VQi(s2) 

<92(si) AQ2(S2) 



Ql{3l)y Ql{32) 
(32(si) AQ2(S2) 



(43) 



and 

^2 



For any number a, 0<77<a<i^<l, we may write a point s € /' as s = + a{s' — s*) 
(s' — s) = (1 — a)(s' — s*). Therefore we can write Q2{s) as Q2{s) = (1 — v)3*{s' ~ s*) +??(! — '7)(s' — 
When s' = 8 we have A = 4 and s« = 4 from Equations (22 23 ) . Both A and s* are strictly increasing as 
functions of s' , so for s' > 8 we have A > 4 and s 



p(si) A I 



> 4. Therefore for s' > 8 we have 
,(S2) = 16(1 -i.2). 



(44) 



It remains to obtain an estimate on the quadratic Qi 



Lemma 6. Assume xq > 9 > 0. Fix the constants < v < rj < 1 and si = + 77A, §2 = s* + j^A with 
A = (s' — s.^,), and s■^, as in Equation (22); let Qi be given by Equation (32) and let s' > 8 V (1 + xf^/9'^). 
Then 

Qi(si) VQi(s2) < A^e^s'. 



Proof. For any number < a < 1, we have a corresponding point s G / for which s = s* + a(s' — s*) and 
(s' — s) = (1 — a)(s' — s*). As shown in Lemma 5j Qi has constant positive second derivative with respect 
to s when s' > 8. Hence its largest value for a E [0, 1] exceeds its largest value for a E [v, 77]. From Equation 
( 32 1 we obtain 



Qi{s)P/ (2^2 



{s' ~ s)ll{s) + s{h{s) - a2) 



= s'l\ — sl\ + sl\ — 2sa2?i + sa\. 
Note that by construction li{s) = hi{s — s*) + 02. Therefore 

Qi(5)/3/(2a2) ^ (5'_5)(5^(5_5^)+a2)2 + ^(5^(^_^^))2 

= [s' - s){hl{s - s^f + 26102(5 - s*) + al) + sbl{s - s^f. 
^For compactness wc henceforth write Qi{s) for Qi{s,s'), leaving the dependence on s' imphcit. 



(45) 
(46) 

(47) 
(48) 
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Write s' — s = (1 — q;)A and s — s* = aA. Then 
Qi{s{a))(i 



2a2 



Qi{s*) 
Qi{s') 



(1 - a)/\{h\a^/^^ + 2&ia2aA + aj) + («' - (1 ^ a)^){cP^^)bl 
A(2(l - a)a&ia2A + (1 - a)al + s'fo^a^A) 
A^((l - a)ae{9^/s' + 1 - .tq) + (1 - a)(e'Vs' + 1 - xq)^ 
+a2As'0V4) 

Qi (s(a)U=o) = 4A(0n/7+1 - xof 
Qi (s(a)|a=i) = A^e^^'. 



(49) 
(50) 
(51) 



For s' > 8 we have A > 4, as shown previously. Therefore for s' > 8 we have the following inequalities 

/A 



7TT < ^< 2 



> 



2^/ZTl 4 



1 > 3- 2^2 > Vs' + 1 - -%/As', 
yZTT - 5\/A7 = 3 - 2^/2 when s' 8 and ^ (V^TT - ^x/A^) = 5^7= - Therefore, 



a;o>6i > dVs' + 1 - ^VA^ 
Qi{s') > Qi(s*). 



(52) 
(53) 
(54) 

(55) 

(56) 
(57) 
(58) 



Inequality (57l follows from (56) because {9\/ s' + I — xq) > 0, by the assumption s' > (l + Therefore 

Qi(si) V Qi(s2) < Qi(s') < A^e^s' 
as was to be shown. □ 



We are now able to state Lemma [7] which leads directly to Theorem [T] 
Lemma 7. Let Xq > 9 > and let s' be sufficiently large that 

»'>«v(l4)v(|!)^ (59) 

Then for the standard Brownian motion given by Equations the conditional first passage time density 

) to arrii 

512 /xo 



g2 defined by Equation (28) for B(s) to arrive at the boundary b(s) given by Equation (16) satisfies 

2 



92 



exp 



32 Vcr 



{s'Y 



(60) 



Proof. Lemmas [s] and [6] and the subsequent remarks establish that when s' > 8 V (1 + xll9'^), the density 
g2 defined by Equation (pSl) satisfies 



92{s') > Ar^{v - rj) exp [-25/(1 - i.^)] 



(61) 
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where the parameters rj and may be chosen arbitrarily subject to the constraints < rj < v < I, and the 
positive terms A and B (which depend on s') are given by 



A{s') 



B{s') 



6Ap0{xo - 0) 
128cr2 



Let us rewrite the right hand side of Equation (61 1 as 172 = AM{ri, ly), introducing 

M{r], v) = T^{v - 77) exp [-25/(1 - v^)] . 



(62) 
(63) 

(64) 



Because the choice of 77 and v is arbitrary, within the constraints, 52 must be bounded below by the supremum 
of M over the region < < < 1. Clearly M is positive and differentiable within this region, and M 
as (77, v) approach any of the boundaries rj = 0, v = 1, -q = v. So the supremum will occur at a point (77+, v+) 
in the interior of the constraint region. Setting dM/dij to zero to find the critical point gives i/+ = 277+. 
Substituting back into ( 64 ) and differentiating gives + B. Substituting again to find 



the maximum value yields 



where we introduce 



Mr, 



1-C 



exp 



-2B 



(65) 



C ^ \/B^ + 2B - B. 



(66) 

Note that for all i? > 0, we have < C < 1. The value of C increases monotonically and C — ^ 1 from below 
as _B — > 00. Since B is directly proportional to s', we can find a value of s' sufficiently large that C exceeds 
any value less than 1. For B > 1/4, for example, we have C > 1/2. Consequently we can bound below the 
exponential factor in (65): 

-IB ' 

> exp [-4B] (67) 



exp 



C 



whenever B > 1/4. 



The factor (1 — C)— S'OasB— >cx). But provided i? > 1 we have the following: 

85^-85 + 1 > 
{4B^+4:B-lf > 16B^{B^ + 2B) 
AB^+iB-l > iBy^ 



2B 



1-C = 1 + B- + 2B > 



AB 



(68) 
(69) 
(70) 

(71) 



Inequality (70) follows from (69) because we clearly have AB + AB — 1 > provided B >1. 



Putting together ([65]), ([67]) and ( [tT] ), we may write Mmax > exp [—AB] / (165) provided B > 1. But requiring 
i? > 1 is equivalent to s' > 128cP/ [(3 A"^ 9^) . As noted earlier, when s' > 8 we have A > 4, so it is enough to 
require that s' > 8cr^/ (/30^) . Therefore under the hypotheses of the lemma, 52(5') > A{s') exp [— 4i?] / (165), 
from which it follows that 



512 fxo \ , 4 

Since A(s') < s' for s' > 0, and since /3 > 0, we have 

512 /xq \ 6 
97(t-1J(^^ 

as required. This completes the proof of Lemma [7] 



/3 /OA 



32 \ a 



I 
32 \ a 



□ 
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3 Strictly positive lower bound for px 

The proof of our main theorem follows from Lemma [7] and Corollary [4] 



Proof of Theorem^ and Remark^ Let {X{t)^t > 0} be an Ornstem-Uhlenbeck process satisfymg Equations 
(T][2 1 and let /3, a and xo be positive constants; let 6* be a constant satisfying xq > 6 > 0. To simplify notation 
we introduce the positive constants 



K 



H = 



512 

l_ 

32 \CT 



- 1 



(72) 
(73) 



Let u = ^ log 



213 




as in Remark 2 Then by virtue of the time rescaling, cf. Equa- 



tion (|10|), from Corollary |4| and LemmalTl provided t > u we may write 

2/3e2^^ 



Px{t) 



> 



92 (e 



2fit 



> 



> 



exp 



1) 

Ij exp 



-Hi 



(74) 
(75) 
(76) 



Let p = 1 + and fc = 2I3K. By hypothesis s = (e^''* - l) > 8 so /3t > log 3 > 1. It follows that 
g-io/3t > exp[-e^^*]. Therefore pjf(i) > 2;9ii:exp [-(1 + iJje^'^*] = A: exp [-pe*^^*] . This concludes the 
proof of both Theorem [T] and Remark [2] □ 



The form of the constants and the bound are consistent with the following heuristic relationships between 
the parameters and the tail of the decay. If /3 is small, the initial condition is far from the threshold and the 
noise is large (xo,(T ^ 0), it is reasonable that the tail of the FPT density should decay relatively slowly. If 
f3 is large, the initial condition is close to the threshold and the noise is small, it is reasonable to expect a 
faster decay of the density at long times. 



4 Discussion 



4.1 Relation to other asymptotic FPT density results 



Nobile et al. [13] derived an asymptotic expression for the first passage time distribution pe{t\xQ) for an 
Ornstein-Uhlenbeck process with limiting mean value Xao to go from initial value xq to a boundary of height 
9 in the limit in which 9 ^ Xq, x^o- They prove that 

pg{t\xQ) = ^e~*^^ + o (7^"*^^) ' as 61 ^ 00, (77) 

where r is the mean first passage time. These same authors extended this result to a broader class of diffusion 
processes in [T3]. Similarly, Giorno et al. showed in [S] that the FPT distribution showed good agreement 
with the form 

p{t) « Z{t)e-^' 
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for passage to a periodically varying boundary 9{t), in cases where 9{t) ^ xq = Xoo- The limit 6 ^ xo,Xoo 
corresponds to the subthreshold case, in contrast to the suprathreshold case xq > 9 > Xoo which concerns us 
here. 



Pauwels [TB] gave conditions on the noise (a) and drift (6) parameters guaranteeing that the stochastic 
differential equation {dXt — <j{Xt)dBt + b{Xt)dt, Xq = xq} should yield a first passage time density, 
p{t, 9\xo), that is jointly continuous and fc-fold differentiable in t, 9 and xq. Pauwels' conditions are satisfied 
by the time homogeneous process considered above. In fS] Lehmann used an integral equation approach to 
extend these results, giving conditions guaranteeing Holder continuity of the FPT distribution. 



Several authors have obtained expressions for moments of the FPT distribution for the time homogeneous 
OU process to constant and some moving barriers, for instance [17| [55] . Lindner [TD] calculated moments 
of the FPT for both exponentially decaying and periodic driving terms similar to the situation described in 
Equation ( 79 ) below. However, knowing the moments of the FPT distribution does not provide a strictly 
positive lower bound for the tail of the density. 



In 1977 Sato proved that the first passage time density pdt) for a Wiener process from B{0) = to a square 
root boundary c^/t + 1 scales as 

p,it) ^ at-P^-^-' (78) 

where < p{c) < 1/2, and it is assumed that c > 0. This provides a lower bound for the tail of the density 
(in the case considered) since if pdt) / (at^^^^-*^^) — > 1 as i — > oo then for any e > 0, there is a T such that for 
t > T we have Pc{t) > at~P^'^'^~^ + e). Sato's result does not provide a lower bound for the suprathreshold 
OUP, however. The suprathreshold case corresponds instead to a boundary of the form — 1 + cy/t after 
appropriate scaling, because in this case the moving boundary crosses the mean value of the Wiener process, 
and the expected first passage time is finite. In the case considered by Sato, the mean first passage time is 
infinite, and the distance from the mean E[B] = to the boundary grows monotonically. 



4.2 Prospective application to neural synchronization 



Strictly positive lower bounds are also of interest in the case of OU processes with time varying forcing. The 
timing of action potentials in nerve cells stimulated by fluctuating current injections are of significant interest 
in neuroscience [6l 1121 126j . There are many models available for theoretical studies of synchronization of 
nerve cells, but for models incorporating the effects of noise the LIF model remains among the most tractable, 
attracting sustained attention [501 1211 121] ■ For this case the model dynamics in Equation [l] is extended to 
incorporate the drive through the time varying function h{t), assumed to be bounded, measurable and of 
zero mean: 

dX = {-/3X + h{t)) dt + a dW. (79) 

When the forcing function h is periodic with period T, one is interested in the succession of first passage times 
{''"ij'''2, ■ ■ • } assuming the process is instantaneously reset to xq upon each encounter with the threshold. It 
can be shown that the sequence of phases of the boundary crossings relative to the periodic drive, (pn = Tn 
mod T, then forms a Markov process on the circle = [0, 1), with transition probabilities given by a map 
K : /fc ^ fk+i where 

/fe+i(0)= / Ki4>,^P)fk{^)d^p. (80) 

Here fk{4>) is the density function for the phase of the fc*'' boundary crossing relative to the periodic drive. 
This framework has been elaborated to study stochastic synchronization of systems analogous to the LIF 
model neuron driven by periodic injected current ^23; as well as LIF model neurons with constant drive 
and periodically varying 9{t) [201 [HI [52] . In both types of models, numerical investigations suggest that 
regardless of initial conditions the phase distribution converges to a unique probability density on the circle 
[20l [23] . This situation contrasts with that of the deterministic periodically forced LIF model neuron, in 
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which a rich collection of resonances and p : q mode locking are known to occur [51 [71 [TSl HT] . The existence 
oip : q mode locked solutions [q boundary crossings per p stimulus periods) implies the existence of multiple 
invariant measures on the circle when p > 1. 



Tateno et al. conjectured that the operator K for the stochastic periodically forced LIF model neuron has 
a unique invariant density and that the sequence of probability densities {K"/o} is asymptotically stable 
f|23j. §5-6). One approach to proving this conjecture would be to exploit the following result (with m — 1): 

Theorem 8 (Corollary 5.7.1 in 8]; see also [53], Equation (38)). Let {X,A,ii) he a measure space, 
if : X X X — M a stochastic kernel, and K the corresponding Markov operator (defined by the integral 
equation {80) above). Denote by Kn the kernel corresponding to K" . If, for some m, 

/ lT^^Kr,^{x,y)dx>{) (81) 

Jx y 

then {-P"} is asympotically stable. 



Exhibiting a strictly positive lower bound analogous to (|3| for the periodically forced case would immediately 
establish the infimum criterion (81 1. The criterion has been assumed to hold under various circumstances (c/. 
Equation (6) in |5; page 330 of [20 ; and Equation (29) in [5T|) but to our knowledge an explicit lower bound 
for the FPT density for a suprathreshold system has not been obtained either for the periodically forced LIF 
model with constant boundary nor for the time homogeneous LIF model with oscillating boundary. The 
extension of Theorem [T] to one or another of these cases remains a topic for future research. 
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